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Abstract 

In the Langevin formalism, the dehcate balance maintained between the fluctuations in the system and their corresponding 
dissipation may be upset by the presence of a secondary, space-dependent stochastic force, particularly in the low friction 
regime. In prior work, the latter was dissipated self-consistently through an additional uniform (mean-field) friction [Shepherd 
and Hernandez, J. Chem. Phys., 115, 2430-2438 (2001).] An alternative approach to ensure that equipartition is satisfied 
relies on the use of a space-dependent friction while ignoring nonlocal correlations. The approach is evaluated with respect to 
its ability to maintain constant temperature for two simple one-dimensional, stochastic potentials of mean force wherein the 
friction can be evaluated explicitly when there is no memory in the barriers. The use of a space-dependent friction is capable 
of providing qualitatively similar results to those obtained previously, but in extreme cases, deviations from equipartition may 
be observed due to the neglect of the memory effects present in the stochastic potentials. 



I. INTRODUCTION 

In the theory of diffusion processes over fixed barri- 
ers, numerous studies have shown that the dissipative 
term in the Langevin equation is rarely constant along 
the reaction coordinateii*2i^*^iS*SiLS A general rate the- 
ory when the friction is both space- and time-dependent 
has been developed to account for this phenomenon over 
the entire friction regime. ^•^''■"'^^ One might naively expect 
that a space-dependent component must be included in 
the friction kernel to capture the essential dynamics of 
a given system. However, this is not always the case. 
Several groups have shown that the average dynamical 
properties may still be adequately described by a gener- 
alized Langevin equation with space- independent friction 
even when the reaction coordinate has a strong spatial 
dependencei24iSiiS An analysis by Haynes and Voth con- 
cluded that the key factor is not whether the friction 
is space-dependent, since it generally will be, but rather 
how the friction varies along the reaction coordinate liS- In 
particular, they suggest that the symmetry of the space- 
dependent friction with respect to the barrier can be used 
as a metric for evaluating the role of the friction in the 
dynamics. Similar product and reactant states will give 
rise to similar (symmetric) friction components about the 
transition state. Perhaps surprisingly, an antisymmetric 
friction does not have a significant impact on the dy- 
namics, while a symmetric friction can result in large 
deviations from the predictions of standard rate theories 
for processes with space-independent frictionA^iSiSiiii^ 
Thus, the Langevin model with a uniform effective fric- 
tion can often approximate the dynamics of projected 
variables even if the formal projection would have re- 
quired a space-dependent model. 

The central question explored in this work is whether 
a single uniform effective friction suffices even when 
the Langevin system is subjected to an external space- 



dependent stochastic potential. The behavior of a 
Brownian particle diffusing across various subsets of 
this class of potentials has been the subject of intense 
research.^^d^d'^d^d^i^°i^b22^23 r^-^is activity has largely 

been motivated by the discovery of resonant activation 
in which the rate of transport over a stochastic barrier 
exhibits a maximum as a function of the correlation time 
in the fluctuations of the barrier height>i^ However, un- 
til recently, simulations of these systems have not been 
performed in the low friction regime, where deviations 
from equipartition may occur, due to an inability to ad- 
equately describe the friction in the presence of an ad- 
ditional stochastic forcc i^^i^^ In previous work, the dis- 
sipation of this excess energy was achieved through a 
self-consistent approach in which the friction constant is 
renormalized iteratively until equipartition is satisfied^ 
This rcnormalization is approximate because it does not 
explicitly account for the correlations between the exter- 
nal stochastic forces across space and time, but rather 
uses a single mean friction to dissipate theses forces at 
times longer than their correlation times. A possible 
improvement to the self-consistent approach can be ob- 
tained by allowing the friction to be space-dependent 
while explicitly ignoring the memory in the stochastic 
potential, In the special case that the stochastic poten- 
tial has no memory, then this treatment is exact. How- 
ever, this approximation is often not justified when mod- 
eling real systems and therefore, the model potentials 
employed are chosen to have an exponentially decaying 
memory of their past states. In the most extreme cases, 
these correlations can result in deviations from equipar- 
tition during the course of the simulation, although the 
space-dependent friction dissipates such fluctuations cor- 
rectly in most situations. The general conclusion appears 
to be that the more detailed space-dependent approach 
is in qualitative agreement with the self-consistent ap- 
proach and hence, as in the fixed barrier case, Langevin 



systems with stochastic forces may be dissipated by a 
single (though renormahzed) uniform friction. 

The conclusions of this work are supported by a study 
of two different classes of one-dimensional problems in 
which the particle diffuses across a periodic array of 
coherent or incoherent barriers. These two cases can 
be specified by sinusoidal or merged-harmonic-oscillator 
potentials, respectively. For such simple forms of the 
stochastic potential, analytic expressions for the friction 
as a function of the spatial coordinate can readily be 
obtained and are presented in Sec. m The resulting 
Langevin dynamics across these potentials dissipated ei- 
ther uniformly or through the space-dependent friction 
are illustrated in Sec. Illll 



II. LANGEVIN MODEL WITH STOCHASTIC 
POTENTIALS 

An equation of motion describing the diffusion of a 
particle influenced by a stochastic potential of mean 
force can be adequately described by a phenomenolog- 
ical Langevin equation of the form. 



and is specified by 



V = -^(t)v + m + F{x;t) , 



(1) 



where F{x;t) = —VxU{x]t) is an external stochastic 
force, and ^(t) is the friction required to dissipate both 
the thermal forces and those due to the external stochas- 
tic potential. The thermal bath is described by ^(i), 
which is a Gaussian white noise source with time correla- 
tion given by the fluctuation-dissipation relation (FDR) , 



mW)) = 2fcBT7th<5(t - t') 



(2) 



In the limit that F{x;t) — F{x;Q) for all t, these equa- 
tions reduce to the Langevin equation with j(t) = 7th- 
Otherwise, the question remains as to what is the appro- 
priate form of j{t). Two approaches for addressing this 
question are presented in Sections |f 1 Bl and lll CI after flrst 
describing the explicit forms of the stochastic potentials. 



A. Stochastic Potential Representation 

The space-dependent friction (SDF) that arises from 
the fluctuations in F(x; t) can readily be evaluated an- 
alytically for two different classes of one-dimensional 
stochastic potentials. The first of these is a sinusoidal 
potential taking the general form, 



U{x;t) 
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(3) 



in which the barriers fluctuate coherently with each other. 
The second is constructed using a series of merged har- 
monic oscillators (MHOs) in which each barrier is allowed 
to fluctuate independently (incoherently) of one another. 
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where the to*'* well and adjacent barrier are centered 



at x^ = —A/2 -I- toA and xj„ = mA, respectively. 
The connection points are chosen to ensure continu- 
ity in the potential and its first derivative such that 
x^ — ±ko\/{2ko — 2fc^) -I- mA. As opposed to the si- 
nusoidal potential, the width of the MHO barriers varies 
stochastically in time according to the relation kj^ — 
— (fco + r]{rn, t)), which, in turn, defines the barrier height 
VX = — fcofc*jA2/(8fco — 8kj^). The remaining parameters 
in the potentials are chosen such that the lattice spacing 
is 4 and the thermal energy of the particle is 1/6 of the 
average value of the barrier heights. 

The stochastic term, ri(t), is defined as an Ornstein- 
Uhlenbeck process governed by the following differential 
equation. 




with the probability distribution, 



and time correlation, 
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The variance of the distribution is given by cr^, Tc is the 
correlation time, and C,{t) is an additional white noise 
source. The distribution of barriers heights for the si- 
nusoidal potential is given directly by the distribution 
of 'q{t), but due to the nature of the expression for the 
barrier heights of the MHOs, the resulting distribution 
for this potential takes on a more complex form that is 
sharper and slightly skewed compared with Eq. El As a 
result, a much smaller range of fluctuations is allowed for 
the MHO than the sinusoidal potential to ensure that the 
distribution does not become significantly non-Gaussian. 
More details on the exact behavior of the MHO barrier 
heights are provided in Ref. l2^ 



B. Uniform Dissipation 

In previous work(24 a self-consistent procedure was de- 
veloped to ensure that the evolution of the system us- 
ing Eq. n remains in thermal equilibrium. This was 
accomplished through an iterative procedure in which 
the friction, given by the sum of the two contributions 
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from the thermal bath and the stochastic potential, i.e. 
7 = 7th + 7f, is renormalized according to the relation, 



v(») 



(8) 



The friction for the next iteration is determined from 
the value of the friction at the current step scaled by 
the magnitude of the deviation from equipartition seen 
in the dynamics until convergence is reached to within a 
desired accuracy. The main criticism to this approach lies 
in the approximation made in developing Eq. |Slin which 
the stochastic potential is treated as a local noise source, 
7f, obeying a fluctuation-dissipation relation equivalent 
to Eq. 121 However, the stochastic potentials have mem- 
ory and are therefore nonlocal in nature leading to non- 
vanishing cumulants at third and higher orders. These 
effects are included, but only in an average manner, to 
second order in this approach. 



C. Space-Dependent Dissipation 

An alternative approach to dissipating the external 
stochastic force relies on replacing the space- and time- 
dependent friction, 7(x, t), by a space-dependent friction, 
7(x(i)), satisfying a local FDR. Given that the size of 
the fluctuations in F(x; t) depend on a; at a given a 
Brownian particle moving quickly across the surface will 
experience a series of forces whose relative magnitudes 
depend on the particle's velocity. However when the the 
Brownian particle moves slowly, the particle will sample 
only the local fluctuations of the stochastic potential in 
the vicinity of its local position x. In this regime, the 
particle arrives at a local quasi-equilibrium which must 
necessarily satisfy the FDR locally. This suggests that 
the dissipation should not be uniform, but rather should 
depend on position, and therefore indirectly on time. It 
should be noted that while the mean-field approach de- 
scribed in the previous subsection is capable of including 
the average of the correlations between the fluctuations, 
the approximation made here does not account for any 
of the memory effects. However, in the limit that there 
is no memory in the external stochastic potential, the 
following results are exact. 

The question now arises of how to explicitly describe 
the friction constant in the presence of an additional fluc- 
tuating force resulting from the potentials of mean force 
given in Eqns. Eland 01 The friction constant must dis- 
sipate the excess energy that arises from the fluctuating 
forces through a local space-dependent FDR, 



(9) 



where the cumulative force is simply the sum of the ther- 
mal Gaussian noise and the stochastic force arising from 
the external potential, Fc — Fth + -Fu- Assuming the 
respective fluctuations in the bath and the potential are 



uncorrelated, i.e. (SFthSFu) 
to 



0, then Eq. ©, reduces 



2fcBT7e(x; t) = m^it')) + {SFvix; tf) , (10) 

The thermal fluctuations are ohmic as given in Eq. 12 
and the relationship for the fluctuations in the force is 
SF\j{x;t) = Fu{x;t) — {Fu{x;t))ij, where the average is 
taken with respect to the auxiliary stochastic variable, 
rj. The average value of the force can be determined 
according to the usual integrals. 
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where the fluctuations in the force are governed by the 
stochastic Ornstein-Uhlenbeck process, ?7, whose proba- 
bility distribution is given by Eq. 

The remaining steps of the derivation rely upon the 
specific form of the potential. As an illustration, the 
SDF is evaluated explicitly below for the simpler sinu- 
soidal (coherent) stochastic potential. (The results for 
the incoherent MHO potential can be found in the Ap- 
pendix.) The derivation begins by direct evaluation of 
Eq. ^Ifor the specific class of potentials. As remarked 
above, the first term reproduces the FDR, Eq. (21 for the 
thermal forces. Ignoring the correlation in the forces at 
different times, the second reduces to: 
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The Gaussian integrals are readily evaluated to yield: 

(5Fu(x;i)2)^4?cos2(^;t). (13) 



Upon substitution into Eq. O the explicit form of the 
SDF is 

2 2 

7c(x;0 =7th'5(i-<') + ^^cos2(^;t) . (14) 

This is the simplest possible form for this result, and is 
due to the separability of the potential into a sum of 
deterministic and linear stochastic terms. In fact, it is 
easily shown that for any separable potential of the form. 



U{x;t) = U{x) +r]{t)W{x) 



(15) 



where U (x) is the deterministic component of the poten- 
tial of mean force, then the additional friction due to the 
stochastic potential is given by 

(SF^ix; tf) = iV,W{x)f f dv {if - 7^)P{ti) , (16) 
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provided the distribution is normalized. The MHO does 
not satisfy the condition of Eq. 1151 and hence its friction 
correction can not be obtained by Eq. 1161 The form of 
the friction correction for the MHO consequently con- 
tains more terms, but the requisite approximation (that 
the forces are uncorrelated at different times) enters the 
derivation in a conceptually equivalent way. 

D. Mean First-Passage Times 

The dynamics of the system were characterized by the 
mean first passage time (MFPT) of a particle to es- 
cape its initial minima and establish a quasi-equilibrium 
within another well. With periodic, stochastic poten- 
tials, this may be accomplished by defining a region of 
the phase space of the particle bounded by an energetic 
constraint. The MFPT is simply the average of a suf- 
ficient number of first passage processes into this region, 
with the corresponding rate given by the inverse of the 
MFPT. While the incorporation of a space-dependent 
friction in the algorithm for the numerical integration 
of the equations of motion would seemingly result in a 
dramatic increase in computational expense, the actual 
effort is comparable to the previous mean-field approach 
because the preliminary convergence procedure for the 
friction constant is now unnecessary. 



III. RESULTS AND DISCUSSION 

The analytic and numerical space-dependent compo- 
nents of the friction over one period of the MHO and 
sinusoidal potentials can be seen in the bottom panel of 
Figs, n and [21 respectively, with the numerical results 
averaged over 500 representative trajectories. 

The top panels display the fluctuations in the poten- 
tial and the resulting forces that give rise to the space- 
dependent friction. The analytic forms of the SDF, dis- 
played as the dotted white line, agree with the corre- 
sponding numerical results, and exact agreement is ob- 
tained upon further averaging. The fluctuations in the 
forces reach a maximum at approximately the midpoint 
between the minima and maxima, where deviations from 
the average force take on the largest values. The fluctu- 
ations in the potential are largest at the barriers, while 
the forces are zero at these locations. This leads to a van- 
ishing contribution to the total friction from the space- 
dependent component at these points. In the well region, 
the behavior of the SDF for the sinusoidal and MHO po- 
tentials is inherently different. The SDF for the MHO is 
zero outside of the barrier region since the wells do not 
fluctuate by construction. However, the sinusoidal poten- 
tial fluctuates continuously throughout leading to a fric- 
tion correction along the entire reaction coordinate. Con- 
sequently, the magnitude of the friction correction in sim- 
ulations employing the sinusoidal potential are slightly 
larger than that in those employing the MHO. But, as 
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FIG. 1: Representative fluctuations over one period of the 
MHO potential and force (top panel), and the resulting space- 
dependent friction (bottom panel). The numerical component 
in the bottom panel is displayed as the solid black line, with 
the analytic result, given in the Appendix, as the dotted white 
line. The temperature is 2/3, the variance is 0.22, and the 
thermal friction is 0.08. 



illustrated below, this effect does not have a dramatic 
effect on the resulting dynamics. 

Values of the friction corrections calculated from the 
iterative and space-dependent approaches for the MHO 
and sinusoidal potentials are displayed in Table ^ with 
the values of the thermal friction listed in the left-most 
column. 

The variance and correlation time for both potentials 
is 0.22 and 1, respectively. The resulting temperatures, 
(fceT = (w^)), are also listed for the space-dependent 
approach. The friction correction in the self-consistent 
method ensures equipartition by definition, and there- 
fore, is not listed. The magnitude of the SDF for all val- 
ues of Tc follow accordingly; however this is the only value 
with respect to the given variance for which any deviation 
from equipartition is observed. As can be seen, both the 
self-consistent and space-dependent components of the 
total friction for each potential provide negligible contri- 
butions for this variance since the magnitude of the fluc- 
tuations in the barrier height are relatively small. There- 
fore the total friction is a sum of a large thermal com- 
ponent, and a space-dependent contribution. The slight 
differences in the magnitudes of the SDF for the two po- 
tentials can be attributed to the piecewise nature of the 
MHO potential. The particles spend most of the simula- 
tion time in the wells which do not fluctuate. A contribu- 
tion to the total friction from the space-dependent term 
is included only when the energetically-limited particle 
accumulates enough energy to explore the upper portion 
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0.08 


0.00 


0.01 


0.67 


0.00 


0.03 


0.69 


0.2 


0.00 


0.01 


0.67 


0.00 


0.03 


0.68 


0.4 


0.00 


0.01 


0.67 


0.01 


0.03 


0.67 



TABLE I: The average of the friction corrections, 7f, calculated by the iterative self-consistent (0) and space-dependent (sdf) 
approaches for the MHO and sinusoidal potentials. The resulting temperatures are also included for the space-dependent 
friction. In all cases the temperature is 2/3 (in units of a standard temperature, fcblo), the variance, — 0.22, and the 
correlation time, Tc = 1. 




FIG. 2: Representative fluctuations over one period of the 
sinusoidal potential and force (top panel), and the resulting 
space-dependent friction (bottom panel). The numerical re- 
sult is displayed as the solid black line, with the analytic re- 
sult, given by Eq. 1141 shown as the dotted white line. The 
parameters used are the same as in Fig. 



of the MHO potential. 

To further explore the accuracy of the space-dependent 
approach, the sinusoidal potential has been studied with 
a ten-fold increase in the variance from 0.22 to 2.2. The 
values of the friction correction from these simulations 
are listed in Table UTI 

The displayed correlation times, Tc, are those that ex- 
hibit the largest resonant activation. Consequently, if 
memory effects in the barrier heights are important in de- 
termining the friction constant, it should be manifested 
here. Although not shown for brevity, outside this re- 
gion of the correlation time, the magnitude of the devia- 
tions from equipartition decrease rapidly, but the size of 
the space-dependent components remains roughly con- 
stant. Similarly, the corresponding corrections arising in 
the self-consistent method also approach zero. As can be 
seen from Table Ull the space-dependent approach results 



in a correction that is roughly constant for all values of 
the correlation time, while the iterative approach does 
exhibit some variation with Tc. This is the expected re- 
sult since the space-dependent friction assumes the fluc- 
tuations in the potential are local and therefore, ignores 
any correlation in the barrier heights. The iterative ap- 
proach, however, is capable of incorporating the memory 
of the potential into the friction correction, but only in an 
average manner. As a consequence, significant deviations 
from equipartition may be observed when simulations are 
performed with a space-dependent friction that ignores 
the correlation effects, as illustrated by this extreme ex- 
ample. 

Figs. 01 and H display the MFPTs obtained for 
the MHO potential with the results from the space- 
dependent and self-consistent approaches in the top and 
bottom panels, respectively. 

The results in Fig. |31 have been calculated using a vari- 
ance of cr^ = 0.05, while those in Fig. ^ use = 0.22. 
The corresponding results for the sinusoidal potential us- 
ing a variance of 0.22 can be seen in Fig. |31 The values 
on the broken axis represent the numerically calculated 
MFPTs in the limits of correlation time, Tc. In the zero- 
correlation time limit, the fluctuations in the potential 
are so rapid that the particle effectively experiences the 
average, stationary potential, from which the dynamics 
were calculated. In the limit of infinite correlations, fluc- 
tuations in the potential are nonexistent, and therefore 
the particle experiences a single realization of the po- 
tential with constant barrier heights determined by the 
initial value sampled from the distribution. The MFPTs 
displayed in Fig. ^obtained with a larger variance alters 
the magnitude of the resonant activation, but influences 
the results for the two approaches equally. The results 
from the simulations with a space-dependent friction are 
systematically shifted to lower MFPTs as seen in all three 
figures. This trend is most readily explained through by 
the trends in Table. ^ In the low friction regime, an 
increase in the friction increases the corresponding rate 
of transport. The average space-dependent contribution 
is always larger than its respective mean field counter- 
part, and is expected to have the largest effect on the re- 
sults with the smallest thermal friction. The fluctuations 
present along the entire reaction coordinate of the sinu- 
soidal potential do not appear to have a dramatic effect 
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0.29 


0.71 
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TABLE II: The average of the friction corrections, 7f, calculated by the iterative self-consistent (0) and space-dependent 
approaches (sdf) for the sinusoidal potential. The resulting temperatures are also included for the space-dependent friction 
method. The temperature is 2/3 in all cases and the variance, = 2.2. 
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FIG. 3: The mean first-passage times (MFPT) for a particle 
diffusing across the MHO stochastic potential are displayed 
for two possible scenarios of the dissipative mechanism. The 
top panel uses space-dependent friction, and the bottom dis- 
plays the uniform friction determined by the self-consistent 
method. The variance for both is 0.05, and the three lines cor- 
respond to values of the thermal friction of 0.08 (solid curve 
with X symbols), 0.2 (dashed curve with triangles), and 0.4 
(dot-dashed curve with squares). The symbols on the broken 
axis represent the numerically calculated MFPTs at the limits 
of the correlation time. 



on the dynamics. The results in Fig. for the sinusoidal 
potential follow the same trend as those in Figs. 01 and 01 
for the MHO potential indicating that the SDF approach 
is capable of adequately describing the fluctuations in the 
system. Aside from the shift, the general behavior of the 
MFPT is adequately reproduced by both methods, par- 
ticularly at larger values of the thermal friction when 
the space-dependent component becomes less significant. 
At this level of description, each of the two approaches 
for constructing the friction are capable of capturing the 
essential dynamics of the system. However, some advan- 



FIG. 4: The mean first-passage times (MFPT) for a particle 
diffusing across the MHO stochastic potential are displayed 
for two possible scenarios of the dissipative mechanism. The 
parameters are the same as in Fig. |21 except the variance is 
0.22. 



tage is gained by using the self-consistent method be- 
cause it ensures the system is kept at constant temper- 
ature for all values of the correlation time throughout 
the simulation, while the space-dependent approach may 
lead to deviations in extreme cases. The most signifi- 
cant difference between the two methods can be seen at 
intermediate correlation times, in which the resonant ac- 
tivation observed from the iterative approach is slightly 
more pronounced. This can particularly be seen in the 
MFPTs when the friction case takes on the smallest value 
of 7th = 0.08. Since the resonant activation arises from 
correlations in the barrier heights, it is not surprising that 
simulations incorporating a friction capable of account- 
ing for this phenomenon can have a noticeable impact 
on the dynamics, even if it does so only in an average 
manner. 
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FIG. 5: The mean first-passage times (MFPT) for a parti- 
cie diflfusing across tlie sinusoidai stociiastic potential are dis- 
played for two possible scenarios of the dissipative mechanism. 
Other than for the change from the MHO to the sinusoidal 
potential, the parameters are the same as in Fig.0] 



IV. CONCLUSIONS 

The space-dependent friction arising from the pres- 
ence of a secondary (external) stochastic potential in the 
Langevin equation has been explicitly derived for two 
simple classes of the stochastic potentials. The numerical 
results are in excellent agreement with analytic expres- 
sions describing the space-dependent friction. The re- 
sulting dynamics have been compared to those obtained 



using an alternate approach in which a uniform correc- 
tion is calculated self-consistently. Although the latter 
approach does effectively include the time correlation be- 
tween the barrier fluctuations at long times, the former 
does not in any sense. This neglect may result in devia- 
tions from cquipartition in some extreme cases. However, 
both approaches are capable of capturing the essential 
dynamics of the system and lead to the now-expected 
resonant activation phenomenon. Consequently, the cen- 
tral result of this paper is that the Langevin dynamics 
of a particle under external stochastic potentials can be 
properly dissipated by a single uniform rcnormalizcd fric- 
tion without loss of qualitative (and often quantitative) 
accuracy. 

The role of the memory time in an external stochas- 
tic potential acting on a particle described by a general- 
ized Langevin equation of motion is still an open ques- 
tion. In this limit, there would presumably be an inter- 
play between the memory time of the thermal friction 
and that of the stochastic potential. When the latter is 
small compared to the former, the quasi-equilibrium con- 
dition central to this work would no longer be satisfied by 
the particle, and hence it is expected that a non-uniform 
(and time-dependent) friction correction would then be 
needed. 
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VI. APPENDIX 



The piecewise nature of the MHO potential results in a piecewise form for the associated SDF. Although incoherent, 
every barrier gives rise to the same averages, and hence the procedure needs to be carried out only over a small region 
defined by the closed interval, [a;5^,a;^j]. The limits of integration over this region can be determined from the 
expression for the connection points 

+m\, (17) 



where kt. 



-{ko 



Iko 2/Crn 

ri(t)). This can equivalently be expressed as 

fcoA 



2ko 



2{xm - mX) 

At the top of the barrier, when x^^ — x}-^, rj{t) = oo. In the intermediate region for arbitrary x, 

fcoA 



(18) 



riit) 



2{x — toA) 



2ko 



V 



(19) 



7 



Otherwise, at the minimum when x,^^ = x^, r]{t) = — fco- 



Although it is apparent from the expression for the barrier height that the corresponding distribution is non- 
Gaussian, the resulting forces are Gaussian with the probability given by Eq. |B| The average force for a given x is 
simply the weighted average of the forces when x is in the respective regions, (x^, a;~ ) and (x~ , x]^), which correspond 
to 7] regions of (— fcoj V*) ^^'^ iv* : The resulting integral for the average value of F{x; t) is now: 



J\dr,F{x)P{r,) + J^drjFix)Pirj) 



(20) 



Here, one must be careful in determining which portion of the force to use in the above equation. For example, when 
77 < rj* , the majority of the force is due to the barrier portion of the potential, not the well component. The average 
can thus be expressed as 



{Fxj{x-t)) 



(21) 



— kr. 



where P'{rf) is defined through the normalization condition, i.e., 



driP'{f]) + / dr]P'{ri) = 1 , 



(22) 



which leads to the probability distribution 



P'iv) 



exp 



(-*) 



1 + erf 



(23) 



where erf(a:) is the standard error function. Use of the normahzation condition reduces the average force to 

{Fv{x;t)) = ko{x - xl) - [ko{x - xl) + ko{x - xl)] / dTjP'{i^) + {x - xl) / dij^jP^Tj) . (24) 
The remaining integrals are readily computed; the explicit form of the average force is 



{Fv{x;t)) = fco(a;-4j - (^o(2a;-4i 



1 + erf 



+ (x~xl,)\j — 





- 


»p(-';cn 


1 + erf 







(25) 



The second quantity to be computed is the average of the square of the force, and the derivation follows that (above) 
of the average force. The limits of integration are the same and the Gaussian integrals can be calculated in the same 
manner. Again using the normalization requirement, the first integral is eliminated such that 



{F^{x-tf) = kl{x-xlf + [kl{x~xlf~kl{x~xlf^ 



+ 2kQ{x - xlf I dr^ijP'iij) + {x- xif 



d-qP'{rf) 



df^-q^P'if^) 



(26) 



— fcn 



The first two integrals are the same as before, and the third can be obtained with little effort. The resulting mean 
squared force is 



(i^u(x;i)') = kl{x~xlf + {kl{x-xlf~kl{x^xlf) 



1 — erf 



1 + erf 



8 



Akoa^ 



t \2 



-I , jx - xj„) 



t ^2 



erf 



exp 



_2 



(-^)-exp(-g) 



erf 



fco 



1 + erf 
77* exp 



(^) 
(-g)+feoexp(-^) 

l + erf(^) 



(27) 



The SDF for the MHO potential is then obtained by appropriate substitutions into Eq. El 
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